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PROGRAM  IN  SHIP  HYDRODYNAMICS 
Executive  Summary  Report 

by 

Robert  F.  Beck,  Project  Director 


The  final  year  of  the  Prognm  in  Ship  Hydrodynamics  has  seen  further  progress  in  the 
understanding  of  the  interaction  of  vorticity  and  the  free  surface.  Budget  cuts  and  a  dels^ 
in  funding  has  caused  reductions  in  the  research  program.  However,  5  Phi>.  theses  and 
one  master  thesis  were  completed  this  year.  The  experimental  research  continued  to 
identify  new  phenomena  using  the  instrumentation  developed  under  this  project  The 
numerical  calculatitms  have  shown  good  i^reement  with  die  experiments  and  fiill  scale 
results. 

Appendix  1  lists  the  year-end  repots  written  the  principal  investigators  for  each  of  the 
research  projects.  These  reports  summarize  die  research  gt^s  and  significant 
accon^^hments  of  the  past  year.  Also  presented  is  pertinent  information  such  as 
publications  and  a  list  of  participants  in  each  project  b  the  secoid  qipendix  are  the 
abstracts  from  papers  published  or  submitted  for  publication  during  the  past  year. 

The  instrumentation  develc^ment  is  continuing.  The  laser  dof^ler  velocimeter  (LDV) 
module  of  the  hydrortynamic  moiitoiing  facility  (HMF)  was  us^  in  experiments  in  die 
towing  tank  for  the  first  time.  The  ID  V  was  mounted  in  a  fixed  position  on  the  side  wall 
of  the  tank  for  the  large  jcxexpeiimoits.  Asubcaniagetocany  dielDVbehindamodelis 
presently  in  the  final  d^ign  phase;  ccmstruction  riiould  start  within  a  mcmth.  Before  die 
LDV  was  mounted  in  the  tank  the  optical  fibers  had  to  be  replaced.  When  the  fibers  were 
pirohased  TSI  had  warned  us  of  th^  finite  life  but  they  did  not  have  enough  mperience 
with  them  to  give  an  accurate  estimate  ci  dm  total  operating  hours  befrxe  replacement  The 
new  fibers  ate  a  ixm-degrading  material  and  shouM  last  in^fmitely . 

The  two  other  subsystems  of  the  HMF  are  also  under  going  modificatitm.  Software  has 
beoi  ^veloped  for  determining  surface  velocities  by  tracking  tagged  fluid  regions.  The 
techniques  ate  similar  to  that  used  in  digital  particle  image  velocimetry.  Originally,  the 
tagged  fluid  regions  were  to  be  warm  spots  created  by  a  scanning  infirared  laser.  Tim 
achievi^le  numba  of  tagged  regions  obtainable  using  ^  technique  is  insufficient  to  get 
the  desired  velocity  and  height  resolution.  Using  particles  to  mmk  the  fiee  surface 
allow  a  much  higher  density  of  velocity  measurements.  To  obtain  height  measurements,  a 
staeo-image  pair  technique  is  presently  being  adapted. 

The  scanning  laser  system  for  multipoint,  refractive  free  surface  slope  nmasurements  is  on 
hold.  Due  to  signal-to-noise  ratio  fnroblems,  the  mepsicron  camera  was  converted  fixnn  a 
fdioton  counting  mode  of  tqmration  to  a  continuous  noode.  During  testing  of  the  new 
system,  the  sensor  tube  failed  Again,  we  had  been  warned  that  the  tube  had  a  finite  life 
but  the  manufacturer  ctwld  not  pre^t  how  Itmg.  The  replacement  cost  is  $26,000  if  a  tube 
can  be  found  because  the  manufacture.  Electro-optics  Division  of  ITT  Corp.,  has 
suspend^  operations.  If  the  stereo-image  system  wcoks,  the  slope  system  will  be 
redundant  It  was  therefore  decided  to  cease  furtimr  developmentuntil  further  information 
is  available. 

A  holognq>hic  imaging  technique  is  under  development  by  Professor  BemaL  The  use  of 
holognq>hy  allows  particles  to  be  tracked  in  three  dimensicms  and  tiius  to  resolve  all  tiiree 
components  of  the  velocity  field  within  a  measurement  volume.  The  technique  has  been 
successfully  used  to  measure  the  velocity  field  sunounding  a  vortex  ring. 
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Research  into  the  fundamental  nature  of  the  interaction  of  vorticity  with  a  free  surface  by 
stu(fying  basic  flows  has  ctmtinued.  Experiments  and  numerical  calculations  widi  tire 
vortex  ring  and  submerged  jet  have  been  carried  out;  particular  emphasis  has  been  place  oa 
the  role  of  surfactants  and  Reynolds  number.  Using  particle  image  velocime^ 
laser  induced  fluorescence  (LI^,  aiul  shadowgraphs  tite  interaction  of  a  vortex  ring  with  a 
contaminated  free  surface  was  examined.  It  was  found  that  the  interaction  varies 
signiflcandy  depencQng  on  the  contamination  number.  For  moderately  to  heavily 
contaminated  surfaces  the  circulation  of  the  secondary  vortices  was  rqrproximately  26%  of 
the  primary  vortex  ring  circulation.  Surface  contamination  also  strongly  influenced  the 
developrrrent  of  three-^Qmension^ty  in  the  vmtex  ring.  Regardless  of  the  omtamination 
numbo;  die  fluid  at  the  surface  was  never  immobilize  Even  for  a  clean  flee  surface,  the 
circulation  of  the  primary  vortex  showol  a  marked  reduction  during  die  interaction  process. 
This  effect  is  believed  to  be  associated  with  vorticity  generated  at  the  surface  by  surface 
curvature  effects. 

Numerical  calculations  f<»^  vortex  ring  interactions  with  the  free  surface  have  been  carried 
out  Fully  nonlinear  free  surface  boundary  crmditions  have  been  used  for  inviscid  flow 
with  a  clean  free  surface.  Non-deformable  or  linearized  flee  surface  boundary  conditions 
were  used  for  the  contaminated  surface  and  viscous  fluid  cases.  Computatitms  agree  very 
well  with  the  experiments.  Both  die  primary  and  secondary  rectmnections  have  been 
simulated.  Full  simulations  using  a  viscous  fluid  and  a  fully  nonlinear  free  surface  are 
presendy  underway. 

Two  sets  of  experiments  have  been  completed  in  the  towing  tank.  The  first  is  the  large 
scale  jet  experiments  using  an  apparatus  Signed  by  Bill  WiUmarth.  The  10  cm  diameter 
jet  is  propeller  driven  and  has  screens  and  a  honeycomb  section  for  flow  straightening  and 
turbulence  management  The  jet  has  a  9:1  area-ratio  exit  contraction.  Measurements 
showed  no  stceamwise  vorticity  exists  at  the  jet  exit  The  large  jet  experiments  are  intended 
to  complement  the  small  jet  experiments  pOTtxmed  last  year  by  Anthony  and  Willmarth. 
Anthony's  results  were  for  a  Reynolds  number  of  12,700  while  the  large  jet  has  a 
nuximum  Reynolds  number  of  300,000.  It  has  been  determined  that  the  late^  spreading 
of  the  surface  current  found  by  Antony  becomes  irulependent  of  Reynolds  numbCT  above 
Re  a  15,000.  Shadowgraph  images  or  the  flee  surface  show  considerable  differerrees  in 
the  flee-surface  signatures  over  the  range  of  Reynolds  and  Froude  numbers  exammed.  At 
high  Reynolds  number  {Re  >100,000)  considerably  more  srrudl-scale  activity  appears  in 
the  free-surface  visualizations.  As  compa^  to  tire  Re  »  12,700  jet,  the  hi^  Reynolds 
number  rrear-surface  layer  is  more  strongly  influenced  by  the  presence  of  the  flw  surface. 

The  second  set  of  experirrrents  in  the  towing  tank  were  continued  tests  with  the  Quapaw 
modeL  Both  propelled  and  unpropelled  cases  have  been  conducted.  As  with  Ae  flee 
surface  jets,  there  is  a  tiun  surface  layer  of  fluid  with  a  strong  outward  (away  flom  the 
wake  centeriirre)  velocity  component  For  the  unpropelled  mottel  stem  vortices  interact  wiA 
the  wake  to  ctny  a  sigmflcant  pmtion  of  the  walre  fluid  downward,  away  from  tire  flee 
surface.  This  interactimi  was  modified  by  the  swirl  of  the  propeller.  Fourier  analysis  of 
shadowgrai>h  images  of  the  wake  behiri^  the  model  has  revealed  a  strong  directional 
component  in  tire  spcetrum  associated  with  the  short  diverging  Kelvin  waves  and  a  broad 
band  spectrum  associated  with  the  waves  at  the  edge  of  the  t^ulent  wake.  A  complete 
arudysis  of  all  the  (Quapaw  results  is  presentiy  being  prepared  for  publication. 

Analytic  and  numerical  nrethods  have  been  used  to  investigate  ntmlinear  waves  in  the  ship 
wake.  Fully  nonlirrear  wave  calculatimis  utilizing  the  desingularized  trrethod  have  been 
employed  u>  investigate  the  solitary  irmer-angle  wave  packets  observed  in  ship  wakes. 
Ctmiparisons  to  lirrear  calculations  using  a  time-domain  Cheen  function  have  revealed  that  it 


is  the  unsteadiness  of  the  flow  field  that  causes  the  wave  packets  and  not  nonlinear 
interactions.  Fully  nonlinear  calculations  have  also  been  done  for  submCTged  bodies. 
Woik  is  proceeding  on  develt^ing  techniques  for  surface  piercing  vessels. 

Analytic  methods  have  also  been  used  to  investigate  the  modification  of  surface  waves  by 
surfactant  slicks.  Surfactant  slicks  cause  surface  tension  gradients  which  in  turn  modify 
the  propagation  of  capillary-gravity  waves.  Results  have  b^  obtained  fcv  the  variation  in 
wave  amplitude  for  &e  two  limiting  cases  of  abrupt  and  slow  spadal  change  in  surfactant 
concentration.  The  methods  developed  are  genoal  so  that  additional  effects  can  be 
incorporated  as  necessary  in  order  to  stu^  realistic  condititms  found  in  the  ocean. 
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APPENDIX  I 


3.1.1.  Ship  Wake  Signatures 

Principal  Investigator:  G.A.  Meadows 

Research  Summary 

Description  of  Scient^te  Research  Goals:  To  understand  the  hydrodynamic 
mechanisms  which  facilitate  the  remote  sensing  of  ship  generated  disturbances, 
experimental  measurements  are  required  which  correlate  the  hydio^manic  prc^)aties  of  the 
flow  field  with  die  electromagnetic  properties  of  the  sensing  &ld.  Efforts  during  t^  most 
recent  period  have  concentrated  upon  more  theoretical  aspects  utilizing  experimentatirm 
conducted  during  previous  portions  of  this  URL  These  efforts  have  hem  dieted  toward 
understanding  the  hydrodynamic  and  electromagnetic  scaling  properties  within  the  surface 
ship  wake. 

The  fundamental  goal  of  diis  research  is  to  investigate  dte  limits  to  u^h  bodi  experimmital 
and  analytical  descriptions  surface  ship  wake  {dienomaia  can  be  utilized  to  determine  the 
hydrocfynamic  sources  of  the  surface  signature.  Theoretical  aspects  of  tins  research  have 
concenuated  <mi  the  solution  of  geophysical  inverse  problems. 

Significant  Results  in  the  Past  Year:  As  a  result  of  previous  experimental 
investigations,  scaling  between  full  scale  and  model  scale  mean  wake  compements  have 
been  well  established  (Meadows  et  al.,  1991).  It  has  been  shown  that  both  the  longitu^nal 
and  transverse  mean  velocity  conqxnients  have  approximately  a  1/x  decay  with  distance 
downstream  based  upon  the  full  set  of  nearly  30,(X]M}  LaGrangian  velocity  observations  in 
the  near  w^  with  correlati<»  coefficients  of  0.98  and  0.94,  respectively.  Ennide  scaling 
was  used  in  the  model  tests.  Preliminary  comparisw  wiA  full  scale  measurements 
indicates  good  agreement  with  the  model  tests  for  both  nagnitude  and  rate  of  decay  of  the 
mean  walm  components. 

Similarly,  the  PIlD.  dissertation  oi  Z.  Wu  has  investigated  scaling  as  used  in  association 
with  the  Kelvin  wake  portion  of  the  ship  wake  signature.  Wu  (1^1)  has  demonstrated  a 
succe^fiil,  constram^  nonlinear  inversion  technique  to  determine  basic  ship  and 
operating  characteristics  from  the  observol  Kelvin  wake  signature.  His  efforts,  again, 
have  utilized  both  laboratory  as  well  as  full  scale  data. 

The  extraction  of  ship  infonnatirm  is  based  on  the  relations  of  the  ship's  wave  spectra, 
wave  anq>litude  fimetion  and  hull  geometry.  In  this  research,  an  analytic  representation  of 
a^ve  elevation  is  introduced  with  the  use  of  tlte  Ifilbort  transform,  and  the  derivation  is 
given  for  tee  calculation  of  the  wave  ampliUute  functirm  from  ^  Fourier  spectrum  of  one 
and  taw  dimensional  complex-valued  wave  elevations.  Methods  and  formulas  are  given 
for  atimating  a  ship's  spe^  and  direction  fix>m  tee  spectrum  of  a  two-dirrtensional  wave 

e*  ,  a  single  wave  cut  oe  two  wave  cuts.  A  theoretical  model  of  tee  wave  amplitude 
on  is  <^elope^  arrd  three  methods  are  designed  for  the  estimatiem  of  a  ship's  length 
firom  ^  wave  arnplitude  functiotL  Under  the  assunq>tion  of  thin-ship  theory,  an  mversion 
technique  to  predict  tee  geometry  of  a  ship's  hull  from,  tee  wave  amplitude  fimetirm  or  its 
magnitude  is  d^eloped  through  the  application  of  a  spectral  meteod  mid  tte  emstrained 
maximum  likelihood  method.  Examples  coiiq)aring  theoretically  data  and  tow 

tank  experimental  data  are  given  to  demonstrate  the  methods  developed  and  estimated 
performance. 

The  final  remaining  component  of  tee  ship  wake  of  interest  is  tee  coitetline  walte.  Our 
previous  investigations  have  indicated  teat  mean  flow  components  decay  rapidly  with 
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distaoce  bdund  tiie  vessel  Hence,  Ae  persisteDce  of  die  late  wake  as  viewed  in  many  SAR 
images  nu^  be  Ae  result  of  large  scale  tuibuloit  vortical  structures  which  persist  for  loi% 
periods  of  time  and  periodiciUy  interact  wiA  the  free  surface.  To  investigate  this 
hypothesis,  Ae  classic  wmic  of  Lamb  (1945)  <»  the  Cauchy-Poisson  disturbance  of  Ae 
surface  was  expan^  u>  three  AmniMons  and  apf^ed  theoretically  to  Ae  surface  sh^ 
wake  probkm.  Preliminary  results  are  encouraging  and  provide  a  plau^le  explanaAm  for 
the  generation  of  small  scale  surface  waves  in  tiie  late  wake.  Results  to  date  of  this  work 
virere  presented  at  Ae  Shq>  Wake  (Zonsortium  workshop. 

LIST  OF  PUBLICATIONS/REPORTS/PRESENTATIONS 

1.  Papers  in  Refereed  Proceedmgs/Ioiimalg 

Wu,  Z.  and  G.  Meadows,  "A  remote  sensmg  technique  for  die  estimation  of  a  moving 
ship's  velocity  and  UmgA  from  its  wave  spectra,”  Proceedings  IEEE  Oceans  '91, 
Hcmolulu,  0(^  1991. 

Meadows,  L.,  G.  Meadows,  A.  Troesch,  S.  Cohen,  KJ*.  Beier,  G.  Root,  O.M. 
GtifBn  and  TJ^.  Swean,  Jr.,  "LaGimigian  Velocity  Profrles  m  Ae  Wake  of  a  I£gh 
Speed  Vessel”  submitted  to  Ocean  Engineering,  1991. 

2.  Technical  Reports.  Non-Refeieed  Papers 

Wu,  Z.,  "On  the  estimation  o(  a  moving  ship's  velocity  and  hull  geome^  information 
from  its  wave  spectra,”  Ph.D.  Thesis,  University  of  Michigan,  Pffl  Technical 
Report  #91-1,  Wl. 

3.  Presentations 

Meadows,  G.,  "Comparison  laboratmy  and  full  scale  wake  measurements,” 
ONR  Code  12  Ship  Wake  Consortium,  1991  Spring  Workshop,  Ann  Arbor,  ^ 
May  1991. 

LIST  OF  HONORS/AWARDS 
none 

LIST  OF  PARTICIPANTS 

O.A.  Meadows,  Associate  Professor 

L.  Meadows,  Research  Associate 
M  Perlin,  Assistant  Prc^essor 

D.  Walker,  Assistant  Professor 
J.  Acuna,  MSE  Candidate 

M.  Song,  PIlD.  m  May  1991 

E.  Wti^fa^  Research  Engineer 
Z.  Wu,  PhJD.  m  August  1991 

D.  Lyzenga,  ERIM 
J.  Lyiden,  ERIM 

R.  Shuchman,  ERIM 

S.  Gaboury,  ERIM 
J.  Schipp,  studmt 
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3.2.6.  Vortical  Structure  of  Ship  Wakes 
Principal  Investigator:  L.P.  Bernal 
Research  Summary 

DtscripHon  of  SeUn^fle  Objectives:  The  main  goal  of  this  project  is  to  obtain  a 
better  amlerstanding  of  the  turbulent  viscous  wake  of  surface  ships  and  underwater 
vehicles  by  means  of  an  experimental  investigation  of  the  underlying  large  scale  vortical 
structure  its  dynamics. 

The  scientific  objective  of  tiie  investigation  conducted  diis  past  year  was  to  obtain  a  better 
quantitative  undmtanding  (tf  the  rede  of  surface  contamination  on  the  (fynamics  vortical 
flows  at  a  £tee  surface.  The  suitability  <df  a  contamination  number  based  on  the  initial 
surface  pimure,  rather  dian  surface  elasticity,  to  scale  the  surface  contaminatitm  effects 
was  investigated. 

A  second  objective  of  tile  research  was  to  develop  Ikdogrqdiic  Particle  Image  Velocimetry 
(HFIV)  for  three-dimeosional  flow  fidd  measurements  in  turbulent  free-8urfiK:e  flows.  The 
main  issue  addressed  in  tiiis  research  was  Ae  resolution  charactnistics  ci  tiie  holographic 
imaging  system. 

Significant  Resuits  in  the  Past  Year:  The  interaction  of  a  vortex  ring  with  a 
contaminated  fiee  surface  was  investigated  to  determine  the  effect  of  surfiKe  contamination 
on  the  d^marnks  of  vortical  flows  at  a  free  surface.  T  jminar  vortex  rings  were  generated 
with  a  piston  type  vortex  ring  generated  propagating  normal  and  toward  tiie  surasce.  Tte 
effect  of  surface  contamination  was  investiga^  ^  depositing  on  the  surface  various 
amounts  of  Oleyl  alcohol.  The  cemtamination  level  was  characterized  by  direct 
measurements  of  surface  tension.  Flow  visualization  of  the  subsurface  flow  using  LDP  md 
the  surface  deformatioo  using  shadowgraph  was  conducted.  The  subsurface  flow  and 
the  fiee  surface  flow  were  mearared  using  PlV.  The  Reimolds  number  range  investigated 

was  ~  »  5,000  -  10,000.  The  contamiiution  number  range  investigated  was 

W  a  a  0  -  310,  where  11 «  surface  pressure;  a  a  diameter  of  tiie  ring;  p  a  coefficiait 

of  viscosity;  F  a  vortex  circulation  strength. 

The  main  ccmclusions  of  this  investigation  can  be  summarized  as  follows:  1)  The 
contamination  number  W  as  defined  ab^  properly  scales  surface  contamination  effects. 
For  W  <  5  tiie  flow  evetiution  is  very  similitr  to  the  evolution  witii  a  clean  free  surface.  For 
W  >  10  the  outward  motion  of  tiie  vortex  ring  is  arrested  by  a  secoiulaty  vratex.  2)  The 
Reynolds  rid|e  is  observed  over  a  limited  range  of  contamination  numbm  5  <  W  <  20. 
For  comparatively  clean  surfaces,  contaminants  did  not  accumulate  in  suffii^t  quantities 
to  form  a  Reynctids  ridge.  For  high  contamination  number,  the  fluid  motion  at  tiie  surfiioe  is 
not  enough  to  clean  tiie  surface  and  form  die  Reyndds  rk^e.  3)  The  circulation  of  the 
secondary  vortices  formed  for  W  >  10  was  measured  using  PIV  and  found  to  be 
qiproxiinuUdy  26%  of  the  initial  vortex  ring  circulation.  4)  The  circulation  (rf  tiie  primary 
vortex  using  PIV  showed  a  marked  reduction  during  tiie  interaction  witii  tiie  surface  even 
fix'  a  ctean  surface.  This  effect  is  believed  to  be  associated  with  vorticity  genoated  at  the 
surface  by  surface  curvature  effects.  5)  Surface  contamination  stnmgly  influenced  the 
developnient  of  three-dimensitmality  in  tiie  vortex  ting.  V<ntex  core  reconnection  to  the 
surface  was  observed  at  low  contamination  numbers.  At  higher  contamination  numbm  the 
formation  of  the  sectmdary  vortex  and  rebounding  of  the  primary  vortex  strongly 
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influenced  the  evolution  of  tiuee-dimrasional  instabilities.  6)  The  fluid  at  die  surface  was 
never  iitimnhiiirerf-  Fiv  measurements  of  die  surface  flow  showed  fluid  motion  even  at  die 
highest  contaminatim  investigated  as  well  as  the  ai^iearance  of  vortical  motions  at  the 
surface. 


In  a  different  investigation,  holographic  imaging  of  small  particles  was  used  to  measure 
three-dimensional  velocity  tields  in  a  volume.  An  in-line  holographic  recording  system  was 
used.  Systematic  tests  were  conducted  to  determine  the  resolution  and  other  imaging 
characteristics  d  hologra(rtuc  systems.  These  experiments  were  conducted  in  water  and  air 
with  particles  of  varicHis  sizes.  A  two  axis  bolo^ldiic  system  was  al^  us^  to  improve 
the  axial  resolution  of  die  measurements.  The  main  conclusions  of  this  investigation  are:  1) 
The  particle  size  resolution  of  a  holographic  system  is  determined  by  die  Fraunhofer 

d2 

number  of  the  holograidiic  recording  system  given  by  F  » where  z  is  the  distance 

fitom  die  particle  to  the  iK^ogrs^hic  plate,  X  is  die  laser  wavelengdi,  and  d  is  the  diameter  of 
particle.  It  was  found  that  the  imaged  particles  must  satisfy  F  >  Fc-  For  holographic 


recording  in  air  without  c^tical  windows  Fc  >  0.003  was  demonstrated.  In  these 
eiqpenments  15  pm  particles  were  consistmtty  imaged  in  a  volume  10x10x10  cnt  The  axial 
position  resolution  ^  a  particle  is  limited  by  die  depth  of  field  of  the  particle  image.  The 
depdi  of  field  was  measured  on  particle  images  and  found  to  be  ±2d^/X.  To  inquove  die 
axml  resolution,  a  two  axis  system  was  deveii^ied  and  used  to  measured  the  velocity  field 
in  a  vortex  ring. 


LIST  OF  PUBLICATIONS/REPORTS/PRESENTATIONS 


1.  Papm  in  Rcfflccd  Precficdings/Jftumalii 

Kachman,  NJ.  and  UP.  Bernal,  "Vortex  ring  dynamics  at  a  contaminated  surface," 
sulnnitted  to  National  Fluid  Psmamics  Congress,  Los  Angeles,  CA,  June  22-25, 
1992. 

Song,  M,  UP.  Beraal,  and  G.  Tryggvason,  "Head  on  collision  of  a  large  vortex  ring 
with  a  fiee  surface,"  tuhaditedvo  Physics  of  Fluids  A. 

Scherer,  J.  and  UP.  Bernal,  "Resolution  characteristics  of  Holographic  Particle  Image 
Velocimetry,"  AIAA  papa  no.  92-0009,  submitted  to  Aeros|^K»  Sciences  Meeting, 
AIAA,  January  1992. 

Kachman,  NJ.,  "The  interaction  d  a  vortex  ring  with  a  contaminated  free  surface," 
PhD  Thesis,  Univenity  <rfh(fichigan,  1991. 

Kachman,  N.J.,  E.  Koshimoto,  and  UP.  Bernal,  "Vortex  ring  interaction  with  a 
contaminated  satfdcc  at  incliiied  incidmice,"  ProceaUngs  Synqxisium  on  Dynamics 
of  Bubbles  and  Vortices  Near  a  Free  Surface,  1991  ASI^  Applied  Mechanics 
Meeting,  Columbus,  June  16-19, 1991. 

Bernal,  L.P.  and  J.  Scherer,  "Holographic  particle  image  velocimetry,"  Bull.  Am. 
Phys.  SoCy  35:10:2238,  1990.  (Abstract  only) 

Kachman,  N.J.  and  UP.  Bernal,  "Interaction  of  a  vortex  ring  with  a  contaminated 
stnface,"  Bull.  Am.  Phys.  Soc,  35:10:2267, 1990.  (Abstract  <»ily) 
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2.  Technical  Rqwrts 


none 

3.  Prcscntatioo& 

Bernal,  LP.  and  J.  Scherer  "3-D  velocity  field  measurement  by  holographic  particle 
image  velocimetiy,"  Advanced  Woficshop  on  Particle  Image,  Speckle  and 
Hologr^hic  Velocimetiy,  Heriot-Watt  University,  Edinburgh  U.K,  Sept  23-25, 
1991. 

Kachman,  NJ.  and  LP.  Bernal,  "Interaction  of  a  vortex  ring  widi  a  contaminated 
surface  at  normal  incidoice,"  ONR  Code  12  Ship  Wake  Conscnrtium,  1991  Spring 
Workshop,  Ann  Arbor,  NO,  May  1991. 

Scherer,  J.  and  LP.  Bernal,  "Recent  experiments  on  holograiriiic  particle  image 
velocimetiy,"  ONR  Code  12  Ship  Wake  Consortium,  1991  Spring  Woikshop,  Am 
Arbor,  MI,  May  1991. 

Bernal,  LP.,  "Experimental  investigations  of  vortex  leconnecticm  at  a  fiee  surface,” 
ONR  1991 Surface  Vmrticity  Woricshop,  UCSD,  La  Jolla,  CA,  Feb.  25-26, 
1991. 

Bernal,  LP.,  "Interaction  of  vortex  rings  witii  clean  and  contaminated  free  surfaces,” 
Invited  Seminar  Department  of  Aerospace  and  Mechanical  Engineering,  University 
of  Arizcma,  Tucson,  AZ,  Feb.  27, 1^1. 

Benud,  LP.,  "A  status  report  on  holographic  particle  in^e  velocimetiy,"  Thermal- 
Fluids  Seminar,  Departnwnt  of  MKhanical  Engineering  and  Applied  Mechanics, 
University  of  Michigan,  Dec.  5, 1990. 

LIST  OF  HONORS/AWARDS 

none 

LIST  OF  PARTICIPANTS 

LP.  Bernal,  Associate  Professor 

N.  Kachman,  PhD.  in  December  1991 

J.  Scherer,  PhD.  Candidate 
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3.3.1.  Nonlinear  Ship  Waves 

Principal  Investigator:  W.W.  Schultz 
Research  Summary 

Description  of  Scientific  Research  Goals:  To  develop  numerical  techniques  that 
will  allow  fully  nonlinear  wave  computations.  Particular  emphasis  is  being  placed  on 
Kelvin  ship  wave  patterns  in  the  far  wake.  Ultimately,  the  model  should  incorporate 
vorticity,  bmaking  waves,  and  subgrid  mottels  of  capilla^  waves  (based  on  the  analytical 
subtasks  of  this  project). 

Significant  Results  in  the  Past  Year:  The  desingularized  boundary  integral 
method  has  been  used  to  study  several  nonlinear  wave  problems  effectively.  The  first 
problem  was  the  generation  of  upstream  runaway  solitons  by  a  moving  disturbance  in 
shallow  water.  We  were  able  to  duplicate  the  results  of  Wu  for  the  weakly  nonlinear  cases. 
We  also  studied  strong  nonlinear  cases  and  the  waves  due  to  submerged  bodies  winch 
caimot  be  considmed  in  Wu's  fKdV  model.  Our  results  were  presented  at  the  Sixth 
Intematicnial  Water  Waves  and  Floating  Bodies  Worksh^.  The  second  problem  was  the 
computations  of  nonlinear  waves  due  to  submerged  bodies.  In  Bertram  et  al.  (1991)  the 
fences  and  moments  on  a  submerged  spheroid  were  computed  using  the  desingularized 
method.  Very  good  agreement  was  found  when  compart  to  both  linear  and  nonlinear 
calculation  me£ods.  The  inner-angle  "solitary"  wavepackets  inside  the  wake  of  a  ship 
were  investigated  using  the  desingularized  method,  the  linear  time-domain  Green  function 
technique  ai^  the  stationary  pha^  aiuilysis.  It  has  been  shown  that  nonlinear  effects  are 
not  essential  for  the  generation  of  the  iimer-angle  wavepackets  while  the  unsteadiness  of  the 
wake  can  result  in  the  innor-angle  wavepackets.  The  abstract  of  a  tentative  paper  entitled 
"Inner-Angle  Wavepackets  in  an  Unsteady  Wake"  has  been  accepted  for  presentation  at  the 
19th  Symposium  on  Naval  Hydrodynamics.  A  manuscript  on  this  topic  is  also  being 
prepared  for  submission  to  Journal  of  Fluid  Mechanics. 

It  has  been  observed  that  surfactant  slicks  cause  appreciable  changes  in  capillary-gravity 
waves.  Widi  this  in  mind,  a  study  has  bemi  carri^  out  by  S.  C^u,  A.F.  Messiter,  and 
W.W.  Schultz  to  predict  the  effect  of  surface  tension  gr^ents  on  the  propagation  of 
Ciq>illaty-gravity  waves.  Analytical  results  for  changes  in  wave  amplitude  have  been 
obtained  in  two  limiting  cases,  for  abrupt  and  slow  spatial  changes  in  surfactant 
concentratimi  and  therefore  in  surface  tension.  For  the  former  case,  solved  using  a 
complex  integration,  a  manuscript  for  publication  has  just  been  complete  For  the  latter 
case,  studied  using  matched  asynq)totic  expansions  and  multiple  scales,  a  manuscript  is  in 
preparation.  While  certain  simpl^ing  assumptions  have  been  used  up  to  this  point,  the 
methods  developed  are  now  avt^able  for  mcorporating  additional  effects  and  fm* 
application  under  realistic  conditions. 

LIST  OF  PUBLICAilONS/REPORTS/PRESENTATIONS 

1.  Papers  in  Rcfaccd  ftwccdinss/Joimials 

Bertram,  V.,  W.W.  Schultz,  Y.  Cao  and  RJ^.  Beck,  "Nonlinear  computations  for 
wave  drag,  lift  and  moment  of  a  submerged  spheroid,"  Ship  Technohgy  Research 
38:3-5,  1991. 

Cao,  Y.,  "Computations  of  nonlinear  gravity  waves  by  a  desingularized  boundary 
integral  mettiod,"  Fh.D.  Thesis,  Univrasity  of  Michigan,  Arm  Arbw,  1991. 
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Joo,  S.W.,  A.F.  Messiter  and  W.W.  Schultz,  "Evolution  of  weakly  nonlinear  water 
waves  in  the  presence  of  viscosity  and  surfactant,"  Journal  ofFlmd  Mechanics, 
229:135-158,  1991. 

Cao,  Y.,  W.W.  Schultz  and  RJF.  Beck,  "Three-dimensional  desingularized  boundary 
integral  methods  fOT  potential  problems,”  Int.  J.  Num.  Meth.  Eng.  11, 1991. 

Schultz,  W.W.  and  J.-M.  Vanden-Broeck,  "Steep  free-suiface  standing  waves,"  Bull. 
Am.  Phys.  Soc,  35:10:2290,  1990.  (Abstract  only) 

Reed,  A.M.,  R.F.  Beck,  O.M.  Griffin  and  RJ>.  Peltzer,  "Hydrodynamics  of  remotely 
sensed  surface  ship  wakes,"  SNAhffi  Transactions  98:319-363, 1990. 

2.  Technical  Reports 

ncme 

3.  Ptesentations 

Cao,  Y.,  W.W.  Schultz,  and  RJ^.  Beck,  "Two-dimensional  solitary  waves  generated 
by  a  moving  disturbance,"  Sixth  Intemadmial  Workshop  on  Water  Waves  and 
Boating  Bodies,  Woodshole,  MA,  1991. 

Messiter,  AB,  "Some  surfactant  effects  cm  deep-water  waves,"  seminars  presented  at 
University  of  East  Anglia,  Ncvwich,  and  at  Inq)erial  College,  University  of  Londcm 
(both  in  England),  Ap^  1991. 

Cao,  Y.,  W.W.  Schultz,  and  RJ’.  Beck,  "Inner-angle  wavq}ackets  in  an  unsteady 
wake,"  accepted  for  presentation  at  19di  Symposium  on  Naval  Hydrodynamics, 
Seoul,  Korea,  Aug.  1992. 

LIST  OF  HONORS/AWARDS 

W.W.  Schultz,  Alexander  vcm  Humboldt  Fellowship,  1991. 

W.W.  Schultz,  University  of  Michigan,  promoted  to  Asscm.  Prof,  with  tenure,  1991. 

LIST  OF  PARTICIPANTS 

W.W.  Scliultz,  Associate  Professor 

R. F.  Beck,  Professcn’ 

AJ^.  Messit»,  Bofessor 

Y.  Cao,  PhJD.  in  August  1991 

S.  Gou,  Bi.D.  Candi^te 
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3.3.2.  Interaction  of  Yorticity  and  Free  Surface  FIo^ 

Principal  Investigator:  G.  Tryggvason 

Rfiscarcli  Summary 

Description  of  Scientific  Goals:  To  contribute  toward  the  understanding  of  the 
interaction  of  vortical  flows  and  the  hree  surface.  Particular  emphasis  is  <m  the  surface 
signature  of  unsteady  flows,  and  how  surface  waves  are  general  and  modified.  This 
understanding  is  sought  by  considering  idealized  model  problems  where  the  various 
mechanisms  may  be  isola^  Generally  the  model  equations,  even  for  fairly  "clean" 
problems,  are  inherendy  n<»ilinear  and  numerical  techniques  are  the  only  feasible  solution 
mediod.  Our  studies  enq)loy  generalized  vortex  methods  and  boundary  integral  techniques 
for  problems  where  viscosity  nu^  be  neglected,  and  fimte  difference  method^  fw  ptoblnns 
whm  it  is  essential  to  account  fw  viscous  effects. 

Significant  Results  in  the  Past  Year:  The  focus  of  our  work  during  last  year  has 
b^  on  the  effect  of  contaminants  <m  viscous  flows  and  wave  generation  due  to  three 
dimensional  flows.  We  have  examined  in  detail  the  collision  of  a  vortex  pair  wiA  a  flat 
(waveless)  free  surface,  fen*  various  amounts  of  surface  cmitaminants.  full  viscous 
fluid  equations  are  stdved  widi  a  finite  differmce  technique.  At  the  surface  the  evolution  of 
the  contaminant  concentration  is  followed,  and  the  sh^  induced  uneven  distribution 
used  as  boundary  conditions  for  the  fluid  motion.  The  results  show  that  relatively  small 
amounts  of  surface  cemtaminant  can  cause  the  vortical  flow  to  behave  as  a  no-slip 
boundary,  in  agreeinent  with  Willmartii's  experiments.  Con:q>utations  fw  various  Reynedds 
numbers  and  equations  of  state  for  the  contaminants  have  been  done.  A  paper  has  been 
submitted  to  the  Physics  Fluids, 

Other  majra:  activities  during  the  last  year  are  presented  in  a  thesis  by  M.  S<»g  (co-chaired 
by  Tryggvason  and  Meadows)  that  was  conq^leted  in  the  spring  of  1991.  Part  of  the  tlwsis 
desml^  experiments  on  a  large  vwtex  ring  in  the  towing  ta^  The  experimental  wo^ 
has  identified  the  nature  of  the  three-dimensional  surface  signatures  that  vppeat  rather  early 
after  the  ring  collides  with  the  surface.  Accompanying  numerical  mortelmg  shows  that  up 
to  the  appeuance  of  three-dimensi<»al  disturbances,  toe  evolution  can  be  explained  by  an 
inviscid  axisymmetric  model.  The  vorticity  is  taken  as  a  single  ring,  and  toe  full  nonlinear 
free  surface  conditions  are  used.  Other  parts  of  the  thesis  focused  on  the  tiiree-dimensional 
interaction  of  both  invisdd  as  well  as  viscous  vOTtex  rings  with  a  free  surface.  Simulattons 
using  an  "essentially"  inviscid  vorttm  model  foUov^  toe  ring  as  it  "opens-up"  (or 
reconnects)  at  die  surface.  Linearized  free  surface  boundary  conditiems  aUow  tte  pcei^tion 
of  the  accompanying  wave  genoation.  In  some  cases  die  simulations  show  a  secondary 
"opening-up"  in  agreement  with  experiments.  Viscous  simulations  using  a  waveless  free 
surface  show  the  same  process,  but  due  to  low  Reynolds  number  only  toe  first 
reconnection  can  be  simulated.  Preliminary  calculations  have  been  used  to  investigate  die 
effect  of  surface  contaminants  on  three-tomensionid  flow.  A  paper  describing  toe 
experiments  and  accompan^^g  modeling  has  been  submitinl  to  the  Posies  of  Fluids^  and 
a  papa  on  die  wave  generation  due  to  vratex  recomwetion,  as  well  as  toe  modding  mme 
con^licated  flows  by  many  "interacti<ms"  where  each  interaction  is  described  by  a  low 
(»der  approximation  a  single  interaction,  is  in  pnqiaratkm. 

Full  simulations  of  the  collision  of  viscous  vratex  rings  with  a  fully  defrxmable  surface  are 
currently  underway.  The  focus  is  on  understanding  toe  vorticity  dynamics  in  terms  of 
baroclinicly  genoated  vorticity  produced  at  the  free  surface.  A  pr^minary  report  was 
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^ven  by  T.  Faical  (self  supported)  at  the  1991  summer  ASME  Apjdied  Mechanics  Meeting 
in  C%io  this  summer. 


A  short  review  of  the  computations  of  vwtex  interactions  with  a  free  surface  and  density 
intmfaces  was  presrated  at  the  AMS/SIAM  summer  seminar  in  Seattle  on  June  18-29, 
1990.  A  manuscript  accompanying  the  talk  will  appear  in  Lectures  in  Applied 
Madiemadcs. 

LIST  OF  PUBLICATIONS/REPORTS/PRESENTATIONS 

1.  Papers  in  Referred  Proceedings/Joumala 

Yu,  D.  and  G.  Tiyggvason,  'The  free  surface  signature  of  unsteady,  two-dimensional 
vortex  flows,"  /.  Fluid  Mech.,  218:  547-572,  1990. 

Tiyggvason,  G.,  J.  Abdollahi-Alibeik,  W.  Willmarth,  and  A.  Hirsa,  "Collisitm  of  a 
vortex  pair  with  a  contaminated  free  surface,"  submitted  to  Phys.  FluidSy  Nov 
1990,  (also  available  as  URI  Technical  Rqwrt  90-1). 

Song,  M.  and  G.  Tiyggvason,  "Free  surface  waves  due  to  the  "opening  up"  of  a  vortex 
ring,"  Bull.  Am.  Phys.  Soc,  35:10:2247,  19M.  (Abstract  o^y) 

Krasny,  R.,  "A  method  for  computing  vortex  sheet  separation,"  43rd  Meeting  of  die 
American  Physical  Society,  Division  of  Fluid  Dynamics,  Nov.  19-21,  Gxnell 
University,  Ithaca,  N.Y.,  Abstracts  in  Bull.  Amer.  Phys.  Soc.  35:2269,  Nov. 
1990. 

Song,  M.,  "Vortex  ring  interacdcms  with  a  free  surface,"  Ph.D.  Thesis,  University  of 
Mich^gaiu  1991. 

Song,  M.,  LJ*.  Bernal,  and  G.  Tiyggvason,  "Head-on  collision  of  a  large  vortex  ring 
with  a  free  surface,"  submitted  to  Phys.  Fluids^  April  1991. 

Faical,  T.,  M.  Song,  S.O.  Unve^,  and  G.  Tiyggvason,  "Collision  of  viscous  vortices 
with  a  free  surfrice  and  dmisity  interfaces,"  Proceedings  of  Symp.  on  Dynamics  of 
Bubbles  and  Vortices  Near  a  Free  Surface",  1991  ASME  Applied  Mechanics 
Meeting,  ed.  Sahin  and  Tiyggvason,  Columbus,  OH,  31-37,  June  16-19, 1991. 

Tiyggvason,  G.,  S.O.  Unverdi,  M.  S<mg,  and  J.  Abdollahi-Alibeik,  "Interaction  of 
vortices  with  a  free  surface  and  demsity  interfaces,"  submitted  to  Lectures  in 
AppB^  Mathematics,  1991. 

2.  Technical  Reports 

none 

3.  BMcmationa 

Tiyggvason,  G.,  "Vortex  Ring  interactiem  with  a  Free  Surface,"  Invited  presentation 
ONR  1991  Free  Surface  Vorticity  Workshop,  San  Diego,  Frt.  25-26, 1991. 

Faical,  T.,  M.  Song,  S.O.  Unverdi,  and  G.  Tiyggvason,  "Collision  of  viscous  vortices 
with  a  free  surface  and  density  interfaces,”  1991  ASME  Applied  Mechanics 
Meeting,  Columbus,  OH,  June  16-19, 1991. 
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LIST  OF  HONORS/AWARDS 


R.  Krasny,  University  of  Michigan,  Faculty  Recognition  Award,  1991. 

O.  ’nyggvascm.  University  of  MicUgan,  im>niolBd  to  Assoc.  Ftof.  with  tenure,  1991. 

LIST  OF  PARTICIPANTS 

Gretar  Tiyggvason,  Associate  Professor 
R.  Krasny,  Associate  Prctfessor 
Musedc  Smg,  FhJ>.  in  May  1991 

lavad  AbdolUihi-Alibeik,  MLS.  Studoit  (no  l<»ger  widt  the  project) 

T.  Faical,  FhD.  Studoit  (self  siq)p(xted) 
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3.3.7.  Turbulence-Free  Surface  biteraction  in  Model  Ship  Wakes 
Principal  bivestigator:  D.T.  Walker 
Research- Symmary 

Description  of  Scientific  Research  Goals:  Tbe  objective  of  this  project  is  to 
determine  the  relati<Miship  between  the  fitee-surface  features  and  the  hydrodynamics  (both 
free-suiface  and  subsurface)  of  model  ship  wakes.  In  the  past  year,  the  scope  of  diis 
project  has  been  increased  to  include  a  m<xe  canonical  flow,  tte  fim-surface  jet  and 
emtinuing  devel<^nnBnt  of  the  HMF. 

Significant  Results  in  the  Pa^  Year:  Progress  has  been  made  on  several  fronts 
since  the  last  report  The  initial  sUufy  of  die  qu^tative  features  of  tbe  wake  of  the  Qutqiaw 
model  have  been  completed.  A  new  large-sc^  jet  qiparatus  has  been  designed  and  tested, 
and  initial  flow  visusdizadon  studies,  as  well  as  turbulent  structure  measurements,  have 
been  conducted.  And  development  of  the  hardware  and  software  has  ctmtinued. 

Our  previous  visualizadcm  studies  of  the  Qu^iaw  wake  have  revealed  die  existence  a  dun 
layer  of  fluid  near  the  free  surface  with  a  strong  outward  (away  from  die  wake  centeiiine) 
vdodty.  This  surface-current  layer  has  been  identified  in  boA  turbulent  surface  jets  ai^ 
wall  jets.  Stem  vortices  which  interact  to  cany  a  significant  portion  of  the  wate  fluid 
downward,  away  from  the  free  surface,  woe  idratifi^  frv  the  unim^lled  model.  This 
interaction  was  modified  by  the  swiri  introduced  by  the  addition  of  a  propeller.  While  the 
detailed  behavior  of  the  stem  vortices  in  the  [uopeUed  case  is  unclear,  visualization  has 
shown  that  these  vortices  no  longer  carry  turbulent  wake  fluid  away  from  the  free  surface. 
This  may  account,  in  part,  for  die  increa^  turbulent  activity  observed  near  the  free  surface 
in  die  wake  of  the  propelled  model.  Shadowgraph  visuaUzatirai  of  the  free  surface  had 
previously  identified  the  existence  of  irregular,  but  somewhat  (uganized,  short  waves  along 
the  edges  of  the  turbulent  wake.  Fourier  an^ysis  of  these  images  has  revealed  a  strong 
directional  component  in  the  spectrum.  Hus  spectral  peak  qipears  to  correspcmd  to  the 
inner  edges  of  tte  steady  Kelvin  wave  system  a^  is  not  related  to  the  waves  at  the  edge  of 
the  turbulent  wake.  A  complete  analysis  of  all  the  Quapaw  results  is  now  being  prt^ued 
for  publication. 

The  turbulent  wake  of  the  (^paw  model  exhibits  many  of  the  same  features  seen  in 
turbulent  jets  near  a  free  surface.  The  jet  has  die  advantage  of  being  experimentalfysinqila* 
to  investigate  than  a  wake,  however  previous  jet  studier  have  been  ctmfined  to  very  low 
Reynolds  number.  To  examine  turbulent  shear  flow  tree  a  free  surface  at  high  Reyrmlds 
number  for  a  mc«e  basic  geometry  than  a  nrodel  shq>  wt.?  u  a  new  large-scale  jet  qqiaratus 
has  been  designed  (by  W.W.  Willnuirth)  and  test^  The  jet  has  a  10  cm  diameter  exit 
nozzle  and  is  capable  cf  operating  at  a  maximum  Reynolds  number  3(X),000.  To  minimize 
cemfinement  effects,  the  jet  is  o^aated  in  the  Naval  Architecture  Department's  large  tow 
tank  (iqiprox.  22  ft  x  10  ft  cross  sectirm).  This  jet  is  propeller  driven  and  has  screens  and 
a  hmeyoranb  section  for  flow  straightening  and  turbulmice  management  and  a  9:1  area-ratio 
exit  contractiem.  Measurements  show  a  maximum  exit  velocity  of  3  m/s  and  a  turbuloice 
intensiQr  of  less  than  0.1  percent  Measurements  showed  no  stieamwise  vratici^  existed  at 
die  jet  exit 

fri  late  June,  a  two-wedc  study  of  the  qualitative  characteristics  of  the  large-scale  jet  was 
undertaken  in  die  large  tow  ta^  For  comparison  purposes,  die  small-scale  jets  used  in  die 
work  of  W.W.  Willmarth  and  D.G.  Anthony  were  examined  at  the  same  time. 
Visualizidiem  of  the  surface  current  showed  that  tire  lateral  spreading  of  the  surfrree  current 
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becomes  independent  of  the  Reynolds  number  above  Re  »  15,000  and  is  independent  of 
Froude  numba.  Conq)aiis<m  with  the  jet  issuing  beneath  a  solid  wall  shows  the  two  types 
of  flow  to  be  qualitatively  similar;  however,  the  surface-current  layer  for  the  free-suiface  jet 
was  thinner  and  spread  more  rapidly  dian  t^  for  the  wall  jet  Shadowgraph  images  of  Ae 
free  surface  show  considerable  differences  in  the  fiee-surface  signatures  over  the  range  of 
Reynolds  and  Froude  numbers  examined.  These  differences,  and  in  particular,  their 
Ftoude-number  dependence,  are  currently  being  examined.  At  die  high  R^olds  numbers 
(>  100,000)  attainable  using  the  new  jet  apparatus,  considerably  mote  small-scale  activity 
appears  in  the  fcee-surface  visualizations.  This  is  due  to  the  larger  range  of  length  and 
velocity  scales  present  at  higher  Reynolds  number. 

Detailed  measurements  of  turbulent  structure  in  the  large-scale  jet  are  have  been  performed 
and  are  now  being  analyzed.  These  experiments  involved  efforts  by  D.G.  Anthcmy  and 
W.W.  Willmarth,  as  well  as  C.-Y.  Chen,  a  graduate  studmt  who  has  recendy  join^  die 
proje^  Measurements  of  turbulence  structure  in  the  large  sc^e  jet  were  acquired  for 
locations  16  and  32  diameters  downstream  of  the  exit  for  a  Reynolds  number  of  101,000. 
These  results  show  considerable  differences  when  compared  with  those  of  Willmarth  amd 
Anthony  in  a  similar  jet  at  Re  -  12,700.  The  near-siuface  layer  is  mote  strongly 
influenoMlby  the  presence  of  die  free  surface  at  the  higher  Reynolds  number.  Thismayfaie 
due  to  the  increased  small-scale  structure  presmit  These  results  are  now  in  preparation  for 
submission  to  the  19^  Symposium  on  Naval  Hydrodynamics  in  Seoul,  Korea. 

For  these  jet  experiments,  the  three-component  laser  velocimeter  developed  by  W.W. 
Willmarth  was  operated  at  a  fixed  location  in  the  tow  tank.  During  testing  of  Ae  laser 
velocimeter  prior  to  installation  in  the  tow  tank,  it  was  found  that  £e  optical  fibms  had 
suffered  sigiuficant  optical  degradation  during  the  last  year.  Degradation  of  the  ff bets  was 
expected,  and  is  a  natural  consequence  of  using  the  instrument  Due  to  the  present  state  of 
knowledge  for  this  type  of  fiber,  the  rate  at  which  this  degradation  would  occur  could  not 
be  accurately  predict^  a  priori.  The  laser  velocimeter  probe  was  returned  to  TSI  Ihc.fQr 
installation  of  new  fibm.  The  new  fibers  are  a  different  non-degrading  material  and 
should  last  indefinitely.  The  new  fibers  are  also  individually  jacketed,  rather  than  bundled. 
This  should  ttiiniiTiize  spurious  velocities  introduced  by  differential  fiber  straining  which 
results  from  vibration  of  the  bundled  fibers.  Design  Im  proceeded  on  a  moving  carriage 
on  which  to  mount  the  laser  velocimeter,  and  construction  wUl  begin  Portly. 

Development  has  proceeded  on  the  two  other  systems  which,  along  with  the  laser 
velocirneter,  comprise  the  HMF.  Software  has  bem  developed  for  determining  surface 
velodties  by  tradong  tagged  fluid  regions.  Due  to  limitations  on  the  achievable  number 
density  of  tagged  regitms  obtainable  using  a  scanning  infirared  laser,  alternative  approaches 
have  been  investiga^  Using  particles  to  mark  the  surface  will  allow  a  Itigher  d^ity  of 
velocity  ineasurements.  The  piwessing  approach  used  determines  the  velocity  field  tom 
pai^cle  displacements  between  two  successive  video  frames.  This  softwiro,  mitially 
^tten  by  L.A.  Meadows,  and  further  developed  by  Ph.D.  student  K.  Thiagarajan, 
inq>lements  an  approach  si^ar  to  that  used  in  ^gital  particle-image  velocimetry.  T^e 
analysis  system  is  currently  being  used  to  examine  subsurface  flow  near  the  edge  of  an 
osdUating  disk.  One  of  the  problems  associated  with  particte-imaging  methods  for  velocity 
measurements  is  getting  consistently  accurate  results.  Spurious  velocity  vectors  appear  in 
the  measured  velocity  field,  seenungly  without  cause.  The  usual  aj^noach  udcen  is  to 
manually  eliminate  these  "bad”  vectors  and  then  fill  tiie  resulting  voi^  in  tlie  data  set  by 
interpolation.  We  are  currently  developing  raticmal  approaches  to  emn'  detection  for  this 
measurement  technique.  Once  these  are  implemented,  the  technique  will  be  tqrplied  to  the 
investigation  of  surface  velocities  to  the  turbulent  jet  fiow  describe  above. 
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This  processing  approach  is  cunently  being  adapted  for  use  in  stereo  image-pair  surface- 
height  measurements.  For  a  stereo  pair,  surface  elevations  are  related  to  the  apparent 
displacement  of  particles  between  the  two  ima^  which  comprise  the  pair.  I]iq)lementati<m 
of  this  technique  requires  accurate  synchron^tion  of  two  video  cameras  mid  real-time 
recording  of  two  simultaneous  images.  These  preliminary  tasks  have  beoi  con^leted  by 
E.A.  Wright  who  has  developed  a  method  fw  stewing  two  numochrome  video  images  in 
one  video  frame  using  color  coding.  The  above  desmbed  analysis  software  for  velocity 
measurements  has  b^  modified  for  this  application  and  testing  of  the  entire  system  is 
underway.  Combination  of  these  two  tqiptoaches  should  yield  hi^-density  measurements 
of  both  surface  velocities  and  elevations. 

Development  is  on  hold  few  the  scanning  laser  system  used  for  multipoint,  refractive  fiee- 
surface  slope  measurements.  The  {wimaty  difficulty  in  making  this  system  both  opo^onal 
and  reliable  has  been  the  imaging  photomultiplier,  or  mepsicron,  u^  for  determining  the 
location  of  the  laser  beam  on  a  screen  suspended  above  the  water  surface.  Due  to  signal-to- 
noise-ratio  problems,  the  mepsicron  was  converted  from  a  photon-counting  mode  of 
operation  to  a  continuous  mode.  This  required  new  processing  electronics  and  data 
acquisition  hardware  and  software.  During  testing  of  these  systems,  the  mqrsicroa  sens<w 
became  inoperable.  At  tUs  time,  it  was  learned  that  the  Electro-optics  Division  of  ITT 
Corp.,  the  manufacturers  of  ^  mepsicrtwi,  had  closed  its  doors.  iW  development,  plus 
tile  mepsicron’s  $26,000  replacement  cost,  have  led  to  a  suspension  of  work  on  this 
system.  If  the  above-outlmed  particle-imaging-based,  stereo-pair  surface  elevation 
measurement  system  proves  successful,  the  slope  measurement  system  will  be  somewhat 
redundant  A  decision  to  resume  work  oa  tte  slope  system  will  require  possibly  the 
adaptation  of  another  sensing  device,  in  place  of  the  mepsicron.  That  decision  will  be  made 
bas^  on  the  outcome  and  evaluation  the  stereo  pair  elevation  measurement  system. 
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Procwdings  IEEE  Oceans  1991 

A  Remote  Sensing  Technique  For  the  Estimation  of  a  Moving  Ship's 
Velocity  and  Length  From  Its  Wave  Spectra 

by 

Z.  Wtt  and  G.  Meadows 

A  lemc^  sensing  technique  to  estimate  a  moving  ship's  speed,  direction  and  length  firom  its 
one  or  two  dim^ional  wave  spectra  is  present^  The  estimati(Mi  of  ship  speed  is  related 
to  the  ship  generated  firee  wave  spectra  and  the  direction  is  related  to  the  positicm  of 
spectrum  loci  in  ^  wave  number  space.  These  relationship  and  the  algorithm  fw  the 
estinution  are  (hscussed.  The  pr^ction  of  ship  length  is  based  on  the  ship  wave 
amplitude  function,  which  can  be  derived  firom  &e  fim  wave  spectra.  Results  of  a 
theoretical  and  experiment  analysis  indicate  that  the  wave  amplitude  function  has  an 
obvious  periodic  character  when  it  is  described  as  a  function  of  the  longitudinal  wave 
num^.  Hence,  ship  length  information,  inherent  in  this  periodicity,  can  recovered. 
The  instantaneous  "nequency"  or  "period"  of  the  wave  amplitude  functicm  may  change 
with  wave  numba  related  to  the  shape  of  the  ship's  bow  and  stem.  Thus,  the  methods  to 
avoid  over-  or  under-estimati<xi  of  ship  Imigth  caused  by  this  frequency  modulati(» 
I^ienomenon  are  also  proposed  together  with  the  analysis  of  the  periodicity.  It  is  shown 
that  high  estimation  accuracy  of  ship  speed  and  lengA  can  be  achieved  for  theoretically 
calculated  wave  data  and  tow  tank  experimental  data.  Relative  errors  are  generally  less  than 
two  percent  for  ship  speed  and  slightly  higher  fat  ship  length. 


In  review.  Ocean  Engineerings  1991 

LaGrangian  Velocity  Profiles  in  the  Wake  of  a  High  Speed  Vessel 

by 

L.  Meadows,  G.  Meadows,  A.  Troesch,  S.  CtAen,  KJ*.  Beier,  G.  Root,  O.M  Griffin, 

and  TJF.  Swean  Jr. 

A  video  imaging  ^hnique  is  enq>loyed  in  a  towing  basin  to  obtain  high  resolution  surface 
LaGrangian  velocitjr  obsavadons  in  the  centoline  wake  region  of  a  high  speed,  twin  screw 
surface  vessel.  Ap^ximately  30,000  surface  velocity  realizations  within  the  wake  region 
extending  from  tte  stem  to  approximately  six  model  ship  lengths  aft  are  used  for 
CQnq>aris<m  wiA  both  recent  full-scale  and  numerical  model  results.  Aiudysis  of  this  data 
set  rev^ds  spadal  decay  rates  in  the  mean  and  fluctuating  velocity  compcments  which  serve 
to  identify  relevant  scate  ci  motion  and  define  candidate  mecbarmms  for  the  persistrace  of 
these  wakes  in  remotely  sensed  ocean  surface  data. 
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Ph.D.  Thesis,  University  of  Michigan,  1991 

On  the  Estimation  of  a  Moving  Ship's  Velocity  and  Hull  Geometry 
Information  from  its  Wave  Spectra 

by 

Z.Wu 

The  wake  generated  by  a  moving  ship  may  extend  fra*  many  tens  of  kilometers  in  the  t^wn 
ocean,  and  can  be  remotely  sens^  Thiou^  indirect  methods,  the  detection  of  a  ship  and 
its  related  characteristics,  is  generally  obtamed  by  measuring  the  ship  generated  waves  or 
their  spec^  From  the  viewpoint  of  remote  sensing,  interesting  problems  exist  related  to 
the  de^tion  of  a  ship’s  presence  and  the  acquisiticm  of  d^amic  and  static  informatirai 
about  it  This  problem  can  be  divided  into  two  basic  aspects.  First  how  to  obtain  a 
moving  ship's  wave  spectra  from  remotely  sensed  images,  and  second,  how  to  extract  dw 
desired  ship  informatum  firan  the  imaged  wave  spectra.  This  thesis  craicentrates  on  the 
latter  aspect  in  particular,  how  to  estimate  a  moving  ship's  direction,  spMd,  Imigdi  and 
hull  shi^  from  its  wave  iqwctra. 

The  extraction  of  ship  information  is  based  rai  the  relations  of  the  ship's  wave  spectra, 
wave  amplitude  function  and  hull  geometry.  In  this  thesis,  an  analytic  representation  of 
wave  elevation  is  introduced  with  the  use  tite  Hilbrat  transform,  and  tlte  derivation  is 
given  fra  the  calculation  of  the  wave  amplitude  function  from  the  Fourier  qpectrum  of  one 
and  two  dimensiraial  coiiq>lex-valued  wave  elevations.  Nfethods  and  formulas  are  given 
for  estimating  a  ship's  spe^  and  direction  from  the  spectrum  of  a  two-dimensional  wave 
patch,  a  single  wave  cut  or  two  wave  cuts.  A  dteoretical  model  of  wave  am|ditude  fuitction 
u  developed,  and  three  meduxls  are  designed  fra  die  estimation  cf  a  ship's  length  from  the 
wave  amplitude  function.  Under  the  assumption  of  thin-ship  the^,  an  inversion 
technique  to  predict  the  geometry  of  a  ship's  hull  ftran  the  wave  anqilitude  function  or  its 
nu^tude  is  developed  through  the  aiqiUcation  of  a  qiectral  method  and  die  constrained 
maximum  likelihood  mediod.  Example  comparing  theoretically  calculated  date  and  tow 
tank  experimental  data  are  given  to  demonstrate  the  methods  developed  and  estimate 
performance. 


Submitted  to  Physics  of  Fluids  A,  1991 
Head-on  Collision  of  a  Large  Vortex  Ring  with  a  Free  Surface 

by 

M.  Song,  LJ*.  Bernal,  and  G.  Ttyggvason 

The  goiraation  of  free  surface  waves  due  to  the  interaction  between  a  free  surface  and  a 
vortex  ring  is  studied  experimentally  and  numerically.  A  well  characterized  vortex  ring 
nramal  to  die  free  surface  is  introduce  The  intraaction  process  is  investigated  by  loddng 
at  the  free  surface  evolution  with  the  aid  of  a  shadowgn^  image,  by  measurement  of  the 
free  surface  elevation,  and  by  visualizing  the  vortex  ring  with  hythogen  bubblc».  Three 
(^erent  free  surface  patterns  are  observed,  deprading  on  the  strmgth  of  the  vortex  ting;  a 
single  circular  depression  of  die  free  surface  fra  weak  vortex  rings  and  very  eariy  times  of 
strong  vortex  r^  interaction,  axisymmetric  radi^  propagating  waves  in  tte  earty  stage  of 
strong  vortex  ring  interaction,  and  frilly  three-dim^sic^  waves  that  depend  on  the  local 
structure  of  d»  vortex  core.  Numerical  simulations  are  used  to  addr^  the  observed 
axrsymmettic  phoiranena. 
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Submitted  to  Aerospace  Scieoces  Meeting,  AIAA  Paper  No.  92-0009 
Resolution  Characteristics  of  Holographic  Particle  Image  Velocimetry 

by 

J.  Scherer  and  UP.  Bernal 

Holographic  Particle  Image  Velocimetry  (HPIV)  uses  holographic  recording  of  a  flow 
seeded  with  small  particles  to  measure  duee-dimcmional  veloci^  flelds.  The  use  of  high 
repetition  rate  lascn  further  extends  the  technique  to  measurement  of  three-dimensimial 
velocity  flelds  as  a  functkm  of  ^e.  The  spa^  resoludoo  characteristics  of  a  typical 
system  are  analyzed.  Basic  requirements  of  die  tologr^diic  recording  and  reccmstru^on 
systems  are  derived.  HPIV  velociQr  measurements  of  a  vrxtex  ring  flow  using  a  two  axis 
holograidiic  recording  system  are  rqxMied  to  illustrate  the  potential  of  the  technique.  Based 
on  these  results  a  spatial  res<duti<m  of  2  mm  and  spatial  dynamic  range  of  100  are  expected 
for  a  typical  HPIV  systent  These  are  comparable  to  typical  resolution  characteristics  of 
LDV  and  (xmventional  PIV  systems. 


PIlD.  Thesis,  Uruversity  of  Michigan,  1991 
The  Interaction  of  a  Vortex  Ring  with  a  Contaminated  Free  Surface 

by 

NJ.  Kachman 

The  interaction  between  a  vortex  ring  and  a  contaminated  free  surface  is  investigated 
expenmentally  using  flow  visualization  and  Particle  Image  Velodmetiy  (PIV).  The  s^ace 
tension  was  measured  directly  on  the  contaminated  surfaces.  A  nondimensional 
contartiination  riumber,  W,  is  int^uced  to  characterize  die  effect  (k  surface  ccntamination. 
hi  die  inclined  incidence  case,  for  low  contamination  numbers,  W  <5  the  flow  evolution 
was  similar  to  that  of  the  clean  free  surface.  AtW  »  9,  the  flow  evolution  changed 
compared  to  die  clean  surface.  In  bottom  vortex  cote  interaction,  tte  vortex  no  longer 
reconnected  to  die  surface.  It  is  shown  that  thiK  was  due  to  the  generation  of  secondary 
yorticity  by  the  surface-active  agents  at  die  surface.  For  all  levels  of  contaminatitxi 
investigated,  the  top  cote  interaction  remained  the  same  as  for  the  clean  case,  widi  die 
vortex  core  reconnecting  to  the  surface.  For  the  normal  incidence  case,  behavior  similar  to 
the  inclined  incidence  case  was  observed.  F^  IV  £  3,  the  evolutitm  of  the  flow  near  the 
surface  was  the  same.  For  10  <  )V  <  20,  the  secondary  vorticity  arrested  the  primary 
vortex  cote  outward  expanskm,  and  forced  it  away  from  tte  surface,  lb  addition  boA  clean 
and  contaminated  rraions  existed  on  the  surface.  For  W  >  20,  the  underwater  flow  was  the 
same  asl0<)V<20butthe  surfue  features  were  gready  re^ed  and  no  cleaned  regions 
on  the  surface  were  produced.  PIV  measurements  are  described  for  the  primary  and 
secondary  vortidty.  A  bias  velocity  was  incorporated  into  the  PIV  to  remove  velocity 
directk)^  ambi^ty.  The  secondary  vorticity  rolled  up  into  a  discrete  vortex  with  a 
circulation  that  increased  for  increasing  W.  The  secondary  vortex  circulatkm  was 
approximately  23%  of  die  primary  vortex.  A  rapid  drop  off  in  circularion  was  noted  to 
occur  whoi  the  ring  reacl^  the  surface.  The  flow  evolution  and  the  amount  of  the 
secondary  circulation  measured  was  similar  to  that  moduced  in  a  solid  wall  interaction. 
PIV  mmauremeats  the  surface  showed  that  the  surrace  velocity  decreased  for  increasing 
contamination  number,  but  for  the  highest  levels  measured  (W  >  360),  die  surface  was 
never  immobilized. 
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Vortex  Ring  Interaction  with  a  Contaminated  Surface  at  Inclined  Incidence 

by 

NJ.  Kachman,  E.  Koshimoto  and  L.P.  Bernal 

The  effect  of  surface  contamination  on  the  interaction  of  a  vortex  ring  with  the  water 
surface  at  inclined  incidence  is  discussed.  A  nondimensional  contamination  number  is 
introduced  to  characterize  the  effect  of  surface  contamination.  At  low  contaminatitm 
number  the  interaction  process  is  very  similar  topologically  to  tiie  interaction  with  a  clean 
surface.  At  larger  values  of  the  contamination  number  tte  topology  of  the  resulting  vortical 
flow  changes.  The  different  behavitv  is  the  result  of  vorticity  generation  at  the  surface  by 
the  surface-active  mateiiaL  It  is  found  that  vortex  reconnection  always  occurs  as  die  upper 
part  of  the  vortex  ring  core  interacts  with  the  free  surface.  In  contrast,  for  contamination 
numbers  higher  than  5,  v<Mtex  reconnection  of  the  lower  part  of  the  vortex  core  does  not 
occur.  The  results  are  also  compared  with  tiie  interaction  of  a  vortex  ring  with  a  solid 
surface. 


Bidl.  Am.  Phys.  Soc.  35:10:2238,  1991 
Holographic  Particle  Image  Velocimetry 
by 

LJ*.  Bernal  and  J.  Scherer 

Holographic  recording  of  small  particles  has  long  been  recognized  as  a  technique  to 
measure  fuUy  three-dimensional  veloci^  flelds.  Ifanne  we  report  the  results  of  several 
ejqwriments  aimed  at  extending  die  technique  to  four-dimenskmal  space-time  velocity  fleld 
measurement  In  the  present  approach  a  Copper  Vapor  Laser  is  us^  to  record  a  sequence 
of  doubly-exposed  holograms  of  the  seeded  flow.  Thus  each  hologram  record  the 
d^lacer^t  of  the  pvticles  at  an  instant  in  time.  The  holograms  are  di^  analyzed  using 
dijptal  image  processing  techniques  to  determine  this  displacement  and  thus  tte  velocity 
fi^d.  Experiments  have  been  conducted  to  determine  fundamental  limitations  of  the 
technique  with  regard  to  ^latial  and  temporal  lestdution  in  Epical  flow  configurations. 


Bull.  Am.  Phys.  Soc.  35:10:2267,  1991 
Interaction  of  a  Vortex  Ring  with  a  Contaminated  Free  Surface 

by 

N  J.  Kachman  and  L.P.  Bernal 

An  experimental  investigation  was  crmducted  to  determine  the  effect  of  different 
concentrations  of  a  surfactant  (oleyl  alcohrd  on  die  Ineaking  and  reconnection  process  of 
vortex  lines  at  a  free  surface.  A  laminar  vortex  ring  with  a  Reynolds  number  of  S(XK)  was 
fired  at  the  fiee  surface  with  an  angle  of  incidence  of  20**.  At  diese  flow  conditions  the 
intoaction  with  a  clean  fiee  sutftee  consists  oi  two  successive  vortex  line  breaking  and 
lecomtection  processes  leading  to  the  formation  of  two  half  vortex  rings  propagating  away 
fiom  each  otfao’.  Flow  visualization  experiments  show  that  the  intoaction  of  die  umer 
vortex  cote  with  the  contaminated  surfitee  results  in  vorticity  generation  and  roll  up.  This 
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does  not  prevent  the  vortex  lines  from  breaking  and  reconnecting  with  the  firee  surface. 
The  interaction  of  the  lower  core  with  the  contaminated  surface  also  results  in  vorticity 
generadcm  and  roll  up  at  the  firee  surface,  bi  dtis  case  however,  the  rolled  up  secondary 
vorticity  causes  the  vortex  core  to  rebound  preventing  vortex  line  breaking  at  the  surface, 
in  addition  a  Remolds  ridge  forms  in  firontof  diis  sec^  interactkm. 


Sh^  Technology  Research  38:3>S,  1991 
Nonlinear  Computations  for  Wave  Drag 
by 

V.  Bertram,  W.W.  Schultz,  Y.  Cao  and  RJ^.  Beck 

The  steady  potential  flow  around  an  elongated  spheroid  under  a  free  water  surface  is 
conq>uted  by  various  linear  and  nonlinear  in^hods.  Results  are  in  good  agreement  and  can 
serve  as  a  validaticm  test  for  future  methods  for  computing  wave  resistance. 


PhD.  Thesis,  University  of  Michigan,  1991 

Computations  of  Nonlinear  Gravity  Waves  by  a  De^ngularized  Boundary 

Integral  Method 

by 

Y.Cao 

A  desingulaiized  boundary  integral  equation  method  ccmbined  wiA  an  Enlerian-T  jgrangian 
time-stepping  technique  is  developed  for  nonlinear  gravity  wave  problems.  The 
desingulanzati<Hi  distance  between  the  boundaty  and  the  sounes  is  related  to  the  local  mesh 
sw  to  ensure  convergence.  Tests  fot  some  simple  problems  show  that  desingularizatimi 
significantly  reduces  the  computer  time  required  to  conq>ute  the  influence  matrix  of  the 
resulting  algebridc  system.  The  algebraic  system  is  still  adequately  well-cmiditioned  to 
allow  fast  iterative  solutions.  Accurate  solutions  can  be  obtained  foe  a  large  range  of 
desingulaiizadim  distances  <m  the  ordo'of  the  mesh  size. 

Several  nonlinear  water  wave  problems  are  dwn  investigated.  The  first  problem  considers 
upstream  runaway  solitons  due  to  a  disturbance  moving  near  critical  speed  in  two- 
dimmsional  shallow  water.  Results  from  die  desingularized  method  with  the  fully 
nonlinear  free  surface  boundary  condition  agree  well  to  those  using  the  fKdV  model  to 
weak  disturbances.  The  fully  nonlinear  model  predicts  larger  sditosthsoi  the  fKdVmodd 
to  strong  disturbances  and  also  predicts  the  breaking  of  waves  to  some  stronger 
disturbances.  Next,  the  problem  of  three-dimensional  waves  due  to  a  submerged  moving 
^herc^  show  good  conqMuison  to  diose  from  other  algorithms.  Finally,  die  generation  of 
irmor-angle  wavepackets  in  the  wake  of  a  ship  is  investigated.  Tl»  tiaoe  most  probable 
causes  of  the  wavepackets  are  examined;  interference  of  &  wave  systems  by  the  bow  and 
stera;  free-surface  nonlinear  effects;  and  wake  unsteadiness  due  to  translation  and 
oscillation  of  the  disturbance.  The  w^  is  studied  with  ntmlinear  calculations  using  die 
desingujarized  mediod  and  with  linear  calculations  using  a  time-domain  Green  function  and 
die  stationary  phase  nwdiod.  It  is  shown  that  nonlin^  effects  are  not  essential  to  die 
generation  and  persistence  of  inner-angle  wavepackets;  die  phenomenon  can  be  exphuned 
by  unsteady  lii^  therxy. 
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Bull  Am,  Phys.  Soc.  35:10:2290,  1991 
Steep  Free-Surface  Standing  Waves 
by 

W.W.  Schultz  and  J.-M  Vanden-Broeck 

A  time-marching  spectral  boundary  integral  method  is  combined  with  a  Newtcm  iteraticm 
procedure  to  se^  spatially  and  temporally  p^odic  solutions  of  standing  waves.  At  lower 
amplitudes,  these  solutions  contibrm  the  existing  perturbation  results.  However,  as  in 
traveling  waves,  we  find  that  the  limiting  form  is  not  the  highest  wave.  We  discuss  the 
limiting  free-surface  accelerations  and  show  results  with  surface  tension. 


SNAME  Transactions  98:319-363, 1990 
Hydrodynamics  of  Remotely  Sensed  Surface  Ship  Wakes 

by 

A.  Reed,  R.  Beck,  O.  GrifGn  and  R.  Peltzer 

Recently,  there  has  been  much  interest  in  the  remote  sensing  of  surface  ship  wakes  using 
s^mthetic  i^ierture  radar  (SAR).  Ship  wake  signatures  obtained  using  SAR  have  several 
<&tinct  features  such  as  a  dark  trailing  centoiine  region,  bright-Une  images  aligned  at  some 
angle  to  the  ship's  path,  and  sometimes  the  I^lvin  wave  pattern.  The  reascms  for  theses 
features  are  not  well  understood.  Many  diffoent  thetnies  are  available  to  explain  <me  ot 
more  of  the  features.  These  theories  are  often  in  conflict  and  no  definitive  answers  are  yet 
available.  Much  analytical  and  numerical  research,  and  both  laboratory  and  full-scale 
experiments  are  presently  underway  in  order  to  develop  jnoper  explanations.  This  paper 
surveys  the  present  research  and  Ae  state  of  the  art  in  the  hydrc^ynamics  of  remotely 
sensed  ship  wakes. 

The  paper  first  describes  in  some  detail  the  distirK^tive  features  of  SAR  images.  Actual 
exarrmles  fiom  early  SEASAT  SAR  images,  and  the  more  recent  JOWIP/Ge<n^  Straight 
and  SIR-B  remote  sensing  experiments  will  be  given.  The  observed  ra^  image 
morphology  includes  dark  centerline  wake  features,  narrow-V  bright  line  features,  and 
Kelvin  wake  manifestations. 

The  discussion  of  the  research  into  surface  ship  wake  hydrodynamics  presently  being 
conducted  at  the  David  Taylm:  Research  Omter,  the  Naval  Research  Laboratory,  and  tiie 
University  of  Michigan  is  divided  into  two  parts:  flow  phenomena  at  and  near  die  free 
surface,  and  submerged  wake  phenomena.  SAR  images  are  formed  by  only  the 
phenomena  right  on  tiie  fiee  surface  since  electrcmiagnetic  waves  do  tKrt  penetrate  the  water 
surface.  However,  these  surface  flows  are  a  numifestation  the  sub-surface  flows  arul  dm 
two  regimes  must  be  investigated  concurtendy.  In  each  part  of  the  paper,  physical 
experiments,  numerical  simulations,  aixi  related  analytical  wmk  are  describe 

Research  into  phenomena  at  and  near  the  free  surface  has  included  investigaaons  into 
surface  films  and  surface-active  materials,  interacti<m  of  turhulence  with  the  fiee  surface, 
and  vortical  flows  including  vortex  ring  and  line  vortex  interactions  with  the  fiee  surface. 
Wave-wave  and  wave-turbulence  interactions,  and  surface  diermal  layos  are  also  topics  of 
active  research.  Experimental,  analytical  and  numoical  studies  have  been  conduct^  and 
some  comparisons  between  the  three  i^^noaches  are  now  possible. 
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While  an  extensive  base  of  submerged  wake  flow  data  is  available  in  the  literature, 
relatively  little  data  has  been  collected  widi  a  free  surface  also  present  This  paper  describes 
recent  research  into  wake  flows  at  and  beneath  a  free  surface,  in  particular  as  it  rehues  to 
surface  ship  wake  hydixxfynamics.  The  ovoall  hydrodynamic  field  which  results  from  die 
ship's  passage  is  discussed  in  terms  of  its  tehdion  to  the  distinctive  surface  signature  which 
can  be  remotely  sensed  by  synthetic  aperture  radar  as  one  principal  example. 


International  Journal  for  Numerical  Metitods  in  Fluids  11, 1990 

Three-Dimensional  Desingularized  Boundary  Integral  Methods  for  Potential 

Problems 

by 

Y.  Cao,  W.W.  Schultz  and  RJ^.  Beck 

The  concept  of  desingulaription  in  three-dimmisional  boundary  integral  computadons  is  re¬ 
examined.  The  boundary  integral  equation  is  desingularized  by  moving  the  singular  points 
away  from  the  boundary  and  outside  die  problem  domain.  We  show  that  the 
desingularization  gives  better  solutions  to  several  problems.  As  a  result  of 
desingularization,  the  surface  integrals  can  be  evaluated  1^  sunpler  techniques,  speeding  op 
the  computation.  The  effects  of  the  desingularization  distance  on  the  solution  and  the 
condition  of  the  resulting  system  of  algebraic  equations  are  studied  for  both  direct  and 
indirect  versions  of  the  boundary  integr^  method.  Computations  show  that  a  broad  range 
of  desin^larization  distances  gives  accurate  solutions  with  significant  savings  in  the 
computation  time.  The  desingularization  distance  must  be  cauefuUy  Imked  to  the  mesh  size 
to  avoid  problems  with  uniqueness  and  iU-amdituming.  As  an  exanq>le,  the  desingularized 
indirect  approach  is  testr^  on  unsteady  nonlinear  three-dimensional  gravity  waves 
generated  by  a  moving  submerged  disturbance;  minimal  computational  difficulties  are 
encountered  at  the  trunoned  boundary. 


Journal  of  Fluid  Mechanics  229: 135-158, 1991 

Evolution  of  Weakly  Nonlinear  Water  Waves  in  the  Presence  of  IHscosity 

and  Surfactant 

by 

S.W.  Joo,  A.F.  Messirer  and  W.W.  Schultz 

A  fonnal  deriyatitm  of  evolution  equations  is  given  for  viscous  gravity  waves  and  viscous 
capillary-gravity  waves  witii  surfactants  in  water  oi  infinite  dq>th.  Multiple  scales  are  used 
to  desoibe  the  slow  modulation  cX  a  wave  packet,  and  noatch^  asymptotic  expansicms  are 
introdiu^  to  represent  the  viscous  boundary  layer  at  the  free  surface.  The  resulting 
dissmative  nonlinear  Schrddinger  equations  show  that  the  largest  terms  in  the  damping 
corfncients  are  unaltered  from  previous  linear  results  up  to  third  order  in  the  anq>litude 
expansions.  The  modulational  instability  of  infinite  wavetrains  of  small  but  finite  amptimrie 
is  studies  numerically.  The  results  show  the  effect  of  viscosity  and  surfactants  on  the 
Benjamin-Feir  instability  and  subsequent  nonlinear  evolution.  In  an  inviscid  limit  for 
capUlary-^vity  waves,  a  small-anq>litude  recurrence  is  observed  that  is  not  directly  related 
to  die  Benjamin-Feir  instability. 
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Journal  of  Fluid  Mechanics  218:547-572, 1990 

The  Free  Surface  Signature  of  Unsteady,  Two-Dimensional 
Vortex  Flows 

by 

D.  Yu  and  G.  Tiyggyason. 

The  inviscid  interaction  of  two-dimensi<»al  vortex  flows  with  a  free  surface  is  studied 
numerically  acombined  vortex/boundary  integral  techniique.  The  vorticity  is  modeled 
as  point  vrutices,  vortex  sheets  and  finite  area  v<»tex  regiras.  Two  problems  are  studied  in 
considerable  dei^,  the  head-on  collision  of  a  v(»iex  pair  with  a  irw  surface  and  the  large 
amplitude  Kelvin-Helmholtz  instability  of  a  submerged  shear-layer.  The  interaction  is 
controlled  by  a  Froude  number  and  by  tte  geomc^  parameters  desoibing  the  initial  vortex 
configuration,  hi  the  large  Frmide  number  limit,  the  surface  motimi  follows  the  vortical 
flow,  but  depends  <mly  weakly  on  the  actual  value  of  the  Froude  number.  Fw  low  Froude 
numbers,  the  free  surface  remains  almost  flat,  and  the  disturbances  caused  by  the  vortical 
flow  decrease  n^idly  with  Froude  number. 


Submitted  to  Physics  of  Fluids  A,  PSH  Tech-  Rpt  90-1, 1990 
Collision  of  a  Vortex  Pair  with  a  Contaminated  Free  Surface 

by 

G.  Tryggvasrm,  J.  Abdollahi-Alibe^  W.W.  Willmarth,  and  A.  Hirsa 

Collisimi  of  a  viscous,  two-dimensional  vwtex  pair  with  a  contaminated,  free  surface  is 
studied  numerically.  The  Froude  number  is  assumed  to  be  small  so  the  surface  remains 
flat  The  full  Navier-Stokes  equations  and  a  ctmservation  equation  for  die  surface 
contaminant  are  solved  numerically  by  a  finite  difference  method.  The  shear  stress  at  the 
free  surface  is  proportional  to  the  ctmtaminaticm  gradient,  arxl  simulati(»s  for  several 
values  of  the  proportionaliQr  constant  (1V)>  ^  as  Reynolds  numbers,  have  been 
performed.  Hw  evolution  is  also  comiM^  with  full-slip  and  no-slip  boundaries.  As  the 
vortices  approach  the  surface,  the  upwelling  between  them  pushes  the  contaminant 
outward,  r^ucing  the  amount  directly  above  the  vortices,  and  lea^g  to  a  clean  region  f<v 
low  TV.  As  IF  is  increased  the  clean  region  bectnnes  smaller,  and  eventually  no  clean 
region  is  formed.  Except  for  very  low  IF,  the  contaminant  layer  leads  to  the  creation  of 
sec(mdary  vortices,  causing  the  miginal  vertices  to  rebound  in  a  similar  way  as  vortices 
colliding  with  a  no-slip  bovmdary.  For  one  case,  the  numerical  results  are  con:q)ared  with 
experimental  measurements  with  satisfactory  results.  Computations  of  a  vmtex  pair 
colliding  obliquely  with  a  contaminated  surface  and  head-on  collision  of  axisymmetric 
vortex  rmgs  are  alro  presented. 
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Bull.  Am.  Phys.  Soc.  35:10:2269,  1991 
A  Method  for  Computing  Vortex  Sheet  Separation 

by 

R.  Krasny 

A  vortex  blob  method  has  been  developed  for  computing  vortex  sheet  separation  and  roll¬ 
up  past  a  solid  boundary  with  sharp  edges.  How  tangency  on  the  bound^  is  enforced  by 
solving  an  integral  equation  for  the  bound  vmrtex  sheet  strength.  A  version  of  the  Kutta 
condition  has  been  implemented  to  determine  the  circulation  shedding  rate  at  an  edge. 
Time-dependent  computadons  are  in  good  agreement  with  Pullin's  results  for  self-simdar 
vortex  sheet  toll-up  (J.  Fluid  Mech.  1978, 88:401).  Several  applications  will  be  presented 
including  die  insU^Uity  of  a  jet  being  expelled  from  a  narrow  slot,  and  the  roU-up  of  the 
wake  behmd  an  oscillating  flat  plate. 


Bull.  Am.  Phys.  Soc.  35:10:2247,  1991 
Free  Surface  Waves  due  to  the  '*Opening-Up'*  of  a  Vortex  Ring 

by 

M.  Scmg  and  G.  Tryggvason 

Experimental  studies  of  the  free  surface  signature  of  underwatm*  vortical  flows  suggest  that 
toe  "opening-up”  of  vortex  filaments  is  a  major  cause  fen*  short  wave  genoadon.  To 
invtttigate  this  phermmena  numerical  simulations  of  toe  collision  of  a  vmtex  ring  of  an 
oblique  angle  of  incident  to  toe  free  surface  have  been  done.  A  nearly  inviscid  flow  is 
assumed,  and  toe  vortex  ring  modeled  by  a  collection  of  vortons.  This  introduces  a  small 
amount  of  dissipation  and  allows  the  ring  to  reconnect  with  its  image  at  toe  free  surface.  In 
most  experi^n^  studies  toe  free  surface  deftnmation  is  very  small,  and  here  toe  free 
surface  rnotion  is  assumed  to  be  linear.  The  geveming  equation  for  the  flee  surface 
deformation  is  solved  in  an  efficient  way  by  a  spectral  method.  Following  the  first 
recoimection  with  toe  surface,  toe  rings  sometimes  recmmect  again,  forming  two  half-rings 
that  propagate  parallel  to  toe  surface.  The  computations  are  compared  to  recent 
experiment^  investigations. 


PhJ7.  Thesis,  University  of  Michigan,  1991 
Vortex  Ring  Interaction  with  a  Free  Surface 

by 

M.  Song 

The  interaction  between  a  vortex  ring  and  a  free  surface  is  studied  experimentally  and 
numerically.  Three  sepmte,  but  relevant  subjects  are  discussed.  First,  a  head-<m  collision 
of  a  large  vortex  ring  with  toe  free  surface  is  investigated  experimentttoy.  Three  diffnent 
free  surface  patterns,  that  depend  on  toe  strength  of  tiie  v(»tex  ring  tmd  the  elapsed  time, 
are  observed:  a  single  circular  depression  of  toe  free  surface  for  weak  vortex  rings, 
axisymmetric  radially  propagating  waves  in  toe  early  stage  of  a  strong  vortex  ring 
interaction,  and  fliUy  thm-dimensional  waves  in  toe  later  stage  fw  both  w^  and  strmig 
rings.  The  early  stage  of  toe  intnaction  is  well  predicted  by  an  inviscid,  axisymmetric 
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numerical  model.  Secondly,  a  generation  of  surface  waves  due  to  a  v<xtex  reconnection  at 
the  free  surface  is  numerically  investigated.  Guided  by  experimental  observations,  Froude 
number  is  assumed  to  be  small  and  hence  free  surface  motion  is  assumed  to  be  linear.  The 
ring  is  modeled  as  a  collection  of  vortons  and  the  free  surface  evolution  is  solved  by  using 
a  spectral  method.  The  proposed  model  successfully  captures  the  reconnection  and  the 
wave  generation,  and  the  results  are  comparable  to  available  experimental  observations. 
The  short  waves  are  generated  as  a  depressed  free  surface  is  released  when  the  vortex  line 
changes  its  topology.  Thirdly,  effects  of  a  surface  contaminant  on  the  oblique  collisicm  of  a 
vortex  ring  with  a  flat  surface  are  studied  numerically.  The  Froude  numba  is  assumed  to 
be  so  sm^  that  that  the  surface  remains  flat  Tbe  full  Navier-Stokes  equations  and  a 
conservation  equation  for  the  contaminant  are  solved  by  a  finite  difference  method.  For  a 
clean  surface,  simulations  at  two  Reynolds  numbers  (200  and  400)  are  done.  For  the 
contaminated  surface,  simulations  with  varying  contamination  parameter  at  the  lower 
Reynolds  number  are  performed.  The  results  show  a  considerable  modification  of  the 
underwater  vortical  flows  due  to  the  surface  vorticity  generated  by  the  contaminant 
gradient  The  distributicm  of  the  ctmtaminant  during  the  interaction  depoids  on  the  strength 
of  the  contamination.  Results  are  also  compared  to  a  collision  with  a  no-sUp  surface. 


Proceedings  Symposium  on  Dynamics  of  Bubbles  and  Vortices  Near  a  Free  Surface 
119:31-37  (Ed.  Sahin  andTryggvason),  ASME  1991 

Collision  of  Viscous  Vortices  with  a  Free  Surface  and  Density  Interfaces 

by 

T.  Faical,  M.  Song,  S.O.  Unverdi,  and  G.Tryggvason 

Full  numerical  simulations  of  the  collision  of  a  viscous  vtxrtex  ring  with  a  fully  deformable 
density  interface  are  presented.  Experiments  and  inviscid  numerical  simulations  of  a  vortex 
ring  or  a  vortex  pair  colliding  with  a  free  surface  or  a  density  interface  show  that  the 
vortices  penetrate  partly  through  the  initial  interface  for  large  Froude  numbers,  but  are 
stopped  when  the  I^ude  numt^  is  small  Fot  small  vortices,  viscous  effects  can  modify 
the  motion,  both  by  reduction  of  the  strength  of  the  v<»tex  ring  itself  and  by  diffusion  of 
baroclirucaily  generated  vrvticity.  For  sni^  Froude  number  Ae  ring  is  deflected  by  the 
interface  and  the  vortex  ring  rectmnects  with  vorticity  created  baroclini^y  at  die  interface. 
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Interaction  of  Vortices  with  a  Free  Surface  and  Density  Interfaces 

by 

G.  Tryggvason,  S.O.  Unverdi,  M.  Song,  and  J.  Abdollahi-Alibeik 

Recmt  investigaticms  of  vortex  interacti(xis  with  a  sharp  density  interface  and  a  free  surface 
^  reviewed.  Particular  emphasis  is  on  numerical  simulations,  but  experimental 
investigations  are  discussed  also.  The  collision  of  a  vortex  pair  and  a  vortex  ring  with  an 
interface  <x  a  free  surface  has  been  studied  in  C(»siderable  detail  In  most  cases  Ae  major 
control  parameter  is  the  Froude  number,  representing  the  relative  strength  of  the  vortices 
and  the  action  of  gravity.  For  high  Froude  numbers  the  vortices  penetrate  through  the 
interface,  and  their  motion  is  only  minimally  nxxiified  for  a  time  that  is  long  relative  to  the 
time  scale  of  the  vortex  motion.  For  low  Froude  numbers,  on  the  other  hand,  the  free 
surface  acts  like  a  rigid  surface  and  inhibits  the  motion  of  the  v(»tices  in  a  direction  normal 
to  the  surface.  Rebounding  is  observed  under  a  variety  of  circumstances,  and  is  found  to 
have  numerous  causes. 
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Observations  of  Turbulence  Near  the  Free  Surface  in  the  Wake  of  a  Model 

Ship 

by 

D.T.  Walker  and  V.G.  Johnston 

Visualization  of  the  flow  in  the  boundary  layer  and  wake  of  a  model  ship  shows  evidence 
of  inboard-rotating  stem  vortices;  however,  the  vortical  structure  persists  only  a  shmt 
distance  aft  of  the  model  before  it  is  effectively  diffused  by  turbulent  mixing.  A  thin, 
rapidly  grovmg,  turbulent  surface  layer,  similar  to  that  seen  in  fiee-surface  turbulent  jets, 
is  also  identifl^  Observations  of  fcw-surface  disturbances  in  the  wake  reveal  that,  near 
the  ship,  there  is  a  region  of  unsteady,  irregular  surface  disturbances  which  correspond  to 
the  subsurface  region  of  turbulent  flow.  At  the  edges  of  tlus  region,  die  small-scale 
disturbances  crea^  by  the  turbulence  coalesce  into  relatively  organizi^  but  unsteady, 
waves.  These  waves  prqragate  outward,  away  fix)m  the  turbulent  region,  and  into  the 
irrotational  free  stream.  Further  aft,  the  distiirbances  caused  by  die  turbidence  become  less 
pronounced;  this  is  due  to  die  gene^  growth  in  die  turbulrat  length  scales  and  decrease  in 
turbulent  velocity  scales  with  distance  aft  At  low  speed,  the  additioic^  a  rotating  propeller 
causes  the  free  surface  in  the  center  portion  of  the  wake  to  exhibit  mark^y  fewer 
disturbances  and  the  swirl  introduced  by  the  propeller  results  in  larger  disturbances  m  <»e 
side  of  the  wake.  At  higher  speed,  the  effects  of  the  propeller  are  much  less  pronounctxl  cm 
the  fiee-surface. 
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Measurements  of  Jet  Turbulence  Beneath  a  Free  Surface 

by 

D.  Anthony  and  W.  Willmarth 

Simultaneous  thiee-con:q)onait  measuiements  d  the  velocity  field  ci  a  submersed  turbulent 
jet  issuing  two  diameters  beneath  and  parallel  to  a  clean  surface  were  made  using  a 
fiber  optic  Laser  Doppler  Velocimeter  ^J3V).  The  measurements  reveal  the  existence  of  a 
strong  outward  current  just  beneath  the  free  surface  that  is  much  wider  than  the  turbulent  jet 
flow  beneath  it  Within  t^  shallow  layer,  turbulent  velocity  fluctuations  rxnrnal  to  die  free 
surface  are  significantly  diminished  while  those  parallel  to  the  surface  are  oihanced.  A 
flow  visualization  experiment  reveals  that  the  surface  current  contains  fluid  structures 
ejected  from  the  jet,  and  that  these  structures  i»opagate  laterally  to  considerable  (flstance 
when  the  water  surface  is  flee  of  surface-active  agmits. 
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Flow  Field  of  a  Submerged  Turbulent  Jet  Interacting  with  a 

Clean  Free  Surface 

by 

W.  Willmarth,  D.  Andiooy  and  A.  Hirsa 

A  0.635  cm  diameter  turbulent  jet  discharging  at  a  sp^  of  200  cm/sec  in  water  parallel  to 
and  2  diameters  beneath  a  clean  fine  surface  is  studied  using  laser  induced  fluorescence 
(UF)  at  Re  »  13,000  and  Froude  number  »  Uexit(gb)*l/2  =  5.66.  Three  component 
turbulence  measurements  show  that  the  jet  spreads  rqiidty  at  a  half  angle  »  35**  with  mean 
velocity  »  10  cm/sec  in  a  thin  layer  near  ^  surface.  LIF  video  images  (to  be  shown) 
reveal  ccdierent  vortical  structures  containing  dyed^  fluid  propagating  outw^  in  this  diin 
layer.  We  si^gest  that  near  the  surface  vertical  mixing  is  inUbited  and  hcuizcmtal  motion  of 
vorticity  parallel  to  the  surface  is  enhanced  owing  to  image  vmticity  above  the  fiee  surface. 
The  rapid  sfueading  near  the  surface  is  suppressed  when  oleyl  alcohol  (a  surfactant  with  9 
dynes/cm  surface  pressure)  is  placed  on  the  surface.  The  suppression  of  horizontal 
spreading  is  attributed  to  secondary  vorticity  generated  when  coherent  vortical  structures 
approach  and  interact  witii  the  surfactant  covert  fiee  surlbce. 
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